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Valley and spin polarization from graphene line defect scattering 
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Quantum transport calculations describing electron scattering off an extended line defect in 
graphene are presented. The calculations include potentials from local magnetic moments recently 
predicted to exist on sites adjacent to the line defect. The transmission probability is derived and 
expressed as a function of valley, spin, and angle of incidence of an electron at the Fermi level being 
scattered. It is shown that the previously predicted valley polarization in a beam of transmitted 
electrons is not significantly influenced by the presence of the magnetic moments. These moments, 
however, do introduce some spin polarization, in addition to the valley polarization, albeit no more 
than about 20%. 
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I. INTRODUCTION 

The success of electronics rests on the ability to con- 
trol electron motion. Such control is typically achieved 
by varying the electrostatic potential in a semiconduc- 
tor with a suitable band gap. This straightforward way 
to control electron motion has proven tremendously suc- 
cessful and is the main reason for the considerable ef- 
fort by the graphene community devoted to graphene 
nanoribbonsir— and bilayer graphene in the presence of 
an electric field4^— Another more subtle way to control 
electron motion is though scattering off deliberate de- 
fects in the material^i^ This approach could add new 
functionality and ultimately prove to be the future for 
electronics. Graphene is a promising material for con- 
trolled electron scattering for several reasons: (i) it has a 
well-defined structure ; 10 ' 11 (ii) owing to its sp 2 hybridiza- 
tion, it has 7r-orbitals near the Fermi level that can form 
extended states^ - (iii) it is a semi- metal with only a 
limited number of scattering channels available near the 
Fermi levelji . and (iv) it offers an electron mobility^ 2 - - — 
high enough to support ballistic transport in the micron 
ranged 

Herein, we consider electron scattering off the extended 
line defect in graphene illustrated in Fig.QJa). This struc- 
ture both preserves the sp 2 hybridization of carbon and 
is precisely defined. It is not a hypothetical structure, 
but one that has already been observed in experiments. 15 
An arbitrary electron scattering off the line defect occu- 
pies a state that away from the line defect approaches 
asymptotically one of graphene, and if the energy of this 
electron is near the graphene Fermi level, it can in addi- 
tion to its energy be identified by its direction of motion, 
its valley, and its spin. Herein, a tight-binding model is 
used to derive the transmission probability of an arbi- 
trary incident electron. The model includes a potential 
to describe ferromagnetically aligned local magnetic mo- 
ments that has been shown to be present in the line defect 
structured These moments break the spin-degeneracy, 
otherwise present, causing there to be spin polarization 



among the electrons of a transmitted beam. This spin 
polarization is found to be rather small, limiting its use. 
More important is that the local magnetic moments do 
not appear to degrade the predicted valley polarization 
of the transmitted electrons near the Fermi levels There- 
fore, the graphene line defect remain an illustrative ex- 
ample of a system where the valley degree of freedom 
can be exploited instead of or alongside the spin degree 
of freedom for applications in quantum information pro- 
cessing. 

The next section develops the theoretical formalism to 
describe the electron scattering off the line defect. This 
theory is applied in Sec. lIIII to the scattering of electrons 
near and at the Fermi level. Conclusions drawn from 
the results are presented in Sec. HV[ including symme- 
try argument explaining the large valley polarization for 
electrons scattering at high angles of incidence. 



II. SCATTERING FORMALISM 

Our objective is to obtain the transmission probabil- 
ity for an arbitrary right-moving electron approaching 
the line defect. The scattering calculations are founded 
on a tight-binding model with a basis set consisting of 
orthonormal 7r-orbitals. 10 Herein, we consider nearest- 
neighbor interactions with a hopping parameter 7 = 
— 2.6eV. Longer-range interactions^ and distortions^ do 
not qualitatively change the results and has therefore 
been ignored for presentational clarity. 

The scattering problem is solved in steps. First, we 
recognize that the structure in Fig.[IJa) can be viewed 
as a set of line defect sites connected to two semi-infinite 
graphene sheets. Next, we solve for the self energy rep- 
resenting all interactions within one of the semi-infinite 
sheets. This can be achieved by exploiting translational 
symmetry along the y-direction. Once the self energy has 
been derived, we can then calculate the retarded Green 
function on the line defect sites, which is needed to obtain 
the sought after transmission probability. 
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FIG. 1. Extended line defect in graphene. (a) An electronic 
Bloch wave approaching the line defect at an angle of inci- 
dence a is being scattered. This scattering is influenced by 
local magnetic moments on the sites (between green arrows) 
next to the line defect sites (between blue arrows), (b) The 
semi-infinite sheet of graphene to the left of the line defect 
transformed into momentum space along the line defect. The 
chain has alternating couplings, 7 and 7', making the self en- 
ergy A representing the influence of the sites away from the 
end point depend on whether the end site is of type A (green) 
or B (blue), (c) Primitive cell of the graphene line defect. 
The sites are paired to form the two chains v = 0, 1. 



To keep the notation tidy, much of the formalism be- 
low is presented in dimensionless units, which can be 
recognized by their diacritic tildes. Energy units are 
scaled with the absolute value of the hopping parameter 
7. Other dimensionless units are also introduced below, 
as needed. 



A. Semi-infinite graphene 

The influence of sites away from the line defect is cap- 
tured by a self energy defined at the edges of the semi- 
infinite sheets of graphene on each side of the line defect. 
To derive this self energy, we first recognize that a semi- 



infinite sheet of graphene with a zigzag edge has trans- 
lational symmetry in the direction along the edge with a 
period equal to the graphene lattice constant a. If a Bloch 
wave is considered with an arbitrary wave vector k y along 
the direction of translational symmetry, we can transform 
the semi-infinite graphene sheet into a semi-infinite lin- 
ear chain with alternating couplings 7 and 7' = 27 cos k y , 
where we have defined k y = k y a/2. See Fig.[ljb). Self- 
energy recurrence relations can be generated through a 
process of replacing the sites away from the end site with 
a self energy, adding a site to the chain, and recalculate 
the self energy at the new end site. In this case, the 
recurrence relations become 



A A = 4 cos 2 k„ IE - A 



A = IE — A 



(1) 
(2) 



where E is the energy of the Bloch wave measured rela- 
tive to the energy at Fermi level and A A is the self energy 
for an end site of type A € {^4, B}. The end site type is re- 
lated to the sublattice of the semi-infinite graphene sheet 
to which the edge sites belong. The retarded solution to 
the self-energy recurrence relations can be expressed as 



A^ = 


1 - 




E - 


A B = 


1 " 




E - 



4 cos 2 k y + 2 cos k y e lkx 



where 



7r + sgn E 



1 + 2 cos k y e 



[ E 2 -1-4 cos 2 L 

/ pi*/ - .. M i. 



(3) 
(4) 



7r — Re arccos 



4 cos k„ 



E 2 - l- 4cos 2 k v 

1m < arccos = - 

I 4 cos k„ 



(5) 



The local moments on the sites adjacent to the line 
defect sites (between the green arrows) can be mod- 
eled through a spin-dependent onsite potential e a , where 
(7 € { — 1,1} is the spin. Assuming the Hubbard pa- 
rameter U = 1.06 [l8| and a difference of the av- 
erage spin population {n a ) per semi-infinite chain of 
(m) — = 1/6, we approximate the onsite potential 

to be e a = U[{n- a ) — (ra CT )] /2 « +0.09 in dimensionless 
energy units for spin a = ±1, respectively^ The local 
moments affect the self energy with the end point at the 
line defect, which is given by 



A 



E-e a - A- 



(6) 



There is translation symmetry, not only in the semi- 
infinite graphene sheets, but also in the full line defect 
structure in Fig. [Ha). The line defect structure has a 
period 2a along the line defect, i.e., twice the period of 
the semi-infinite graphene sheet. The primitive cell of 
the line defect structure is shown in Fig.[TJc), where the 
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sites have been divided into bottom, v — 0, and top, 
v = 1, sites. Rather than using the real space basis \v), 
it is more convenient to use a basis \n), in which the self 
energy is diagonal. To find this basis, we first recognize 
that the wave vector of the line defect structure must 
be conserved in the scattering process. As a result, the 
Bloch wave in the semi-infinite graphene sheet with wave 
vector k y can only couple to one other Bloch wave in the 
scattering process, the one with wave vector k y + n/a. 
From the phase relationship between equivalent sites with 
u — 0,1 imposed by the translational symmetry of the 
semi-infinite graphene sheet, we obtain 



i»> = 4=E 



= i(2fcy+rt7r)y 



(7) 



where n = 0, 1. As the states |n) are eigenstates of the 
self-energy operator S, we have 



(nlEln') = £„<W, 



(8) 



where E n is the self energy A in Eq. calculated for 
the Bloch wave with wave vector k y + nit / a. 

In calculating the transmission probability, we also 
need the elements of the broadening operator T = 

i ^E — . From this definition, we sec that \n) are 
also eigenstates of T, yielding the elements 



(n\f\n') = -2 Im 



(9) 



Line Defect 



Let us focus on the center two sites in the primitive cell 
in Fig. [3c) forming the line defect. With the coupling of 
these sites to all other sites in the primitive cell already 
accounted for through the self energy, the only coupling 
in the Hamiltonian is the coupling 7 between the two 
sites with v = and v = 1. In the basis set defined by 
Eq. (J7]) , the Hamiltonian elements are 



(n\H\ri) 



gi(2k y +n'n) _|_ g — i(2k y +n-ir) 



(10) 



Given the self energy and the Hamiltonian, we can cal- 
culate the retarded Green function operator G = (EI — 

H — 2S) , where / is the unit operator. The factor 2 
in front of the self energy operator reflects the connec- 
tions of the line defect to the two semi-infinite graphene 
sheets. Using Eq. © and Eq. <jTUJ) , we find the elements 
of the retarded Green function operator, which can be 
expressed as 



(n\G\ri) = < 



B-cos(2fe y +n7r)-2Si_ r , 
det (itl-H-2^) 

i sin(2fc y +ro7r) 
dot (EI-H-2t) 



for 11 = 



for li/li', 



(11) 



where 

det (EI 



H - 2S) = (E - 2E ) (E - 2Ei) - 1 

fh). (12) 



2cos2A: y (Eo 



Using the elements of the broadening operator in Eq. (|9]) 
and the Green function operator in Eq. (fl"2|) , we find that 
the probability that a state \n) transmits through the line 
defect into state |n') is given by 



T n ^ n , = (n|f \n)(n\G\n')(n'\f\n')(ri\&\n) 



4 Im XL Im S r 



[n\G\n') 



(13) 



III. ELECTRONS AT THE FERMI LEVEL 

Because the scattering process conserves energy and 
the wave vector associated with the translation symme- 
try along the line defect, it is natural to develop the scat- 
tering formalism based on these parameters. Although 
these parameters together with the index n form a pa- 
rameter space covering all possible asymptotic graphene 
states, there are other more intuitive representations. 
An arbitrary graphene state is typically described by a 
spin, a band index, and a two-dimensional wave vector 
k = (k x ,ky). Rather than using the conventional hexag- 
onal first Brillouin zone in graphene, it is herein more 
convenient to use an alternative reciprocal primitive cell 
bounded by k x £ [0, 4-7r/-\/3a[ and k y € [— 7r/a, 7r/a[. As- 
suming an extended right-moving asymptotic graphene 
state, the wave vector component k x is determined from 
Eq. (|5|) with k x = \f%k x aj2. The graphene band struc- 
ture, given byi2> 



E = rj a / 1 + 4 cos 2 k y + 4 cos k y cos k x , 



(14) 



where r\ = ±1 refers to the conduction and valence bands, 
is shown in Fig.[2]Ja). To make the asymptotic graphene 
state commensurate with the line defect structure, we 
fold the graphene band structure as shown in Fig.[5Jb), 
where in this case k v € {—ir/2a,ir/2a[. The bands fur- 
thest and closest to the Fermi level, located where the 
Dirac cones meet, correspond to n = and n = 1, re- 
spectively. 

Let us focus our attention to the asymptotic states 
of most interest, i.e. those near the Fermi level. For 
energies with \E\ < \J~2~ — 1, which include all energies 
within an eV from the Fermi level, ImXo = 0, i.e. all 
extended states are in the n = 1 band. In this energy 
regime, elastic scattering does not permit any interband 
scattering at the line defect, and thus the transmission 
probability T T a , given by T\^i in Eq. (fT3"|) , is 



l(lm£i) 2 (E + cos 2k y - 2S " ^ 



(E - 2S ) (E - 2Si) - 1 + 2 cos 2k y (S - £1 

(15) 



(b) 




FIG. 2. Energy of the asymptotic states given by the graphene 
band structure, (a) The band structure presented using a 
rectangular reciprocal primitive cell, in which the the two val- 
ley t = ±1 are shown as blue and red, respectively. The en- 
ergy is set to zero at the Fermi level, (b) Graphene band struc- 
ture folded to be commensurate with the line defect structure. 
Bands farthest and closest to the Fermi level have n = and 
n = l, respectively. 



The expression above is an exact result that can be eval- 
uated numerically. To make further analytical progress, 
we focus on the low energy limit. In doing so, we first 
introduce the graphene wave vector q = (q x ,q y ), where 
1x = k x — 2ir /v3a and q y = k y + 2nr/3a, centered at 
valley r G { — 1,1}. In our low-energy regime, we can 
use r, a, and q to describe the asymptotic graphene 
state. To lowest order in q = \q\, the energy disper- 
sion in Eq. (fT4"|) is E = rjhvpq, where vp — y/3\j\a/2h is 
the Fermi velocity. Next, we introduce the angle of inci- 
dence a shown in Fig.[IJa). From the group velocity rela- 
tion tana = (dE/dq y ) j (dE/dq x ) = q y /q x and the as- 
sumption of a right-moving state, which gives sgn q x = tj, 
we obtain q x = (E/Hvf) cos a and q y = (E/Hvp) sin a. 
Therefore, the asymptotic graphene state can be ex- 
pressed uniquely by the energy E, valley r, spin er, and 
angle of incidence a of the incident electron. After ex- 
pressing the transmission probability in Eq. (|15[) using 
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FIG. 3. Transmission probability of an incident electron at 
the Fermi level in valley r and with spin a, approaching the 
line defect at the angle of incidence a. 



these quantities, we find the zero energy limit 

2 

Im Si 



1 + Si 



1 



1 



1— sin ra+Scr (l— 2 sin ra j 



From this expression, shown in Fig. [31 we see that 
T TjfT (a) = T a (ra), which implies that changing valley 
has the same effect as changing the sign of the angle of 
incidence. By letting dT^/dra — 0, we find that the 
transmission has a stationary point at 



ra 



l + 2e CT 



1 



(17) 



Inserted into Eq. (|16j) . this relation gives the maximum 
transmission 



fpmax 



_ / (1 - 2e ff - e 2 a + 2e 3 a + e 4 a ) 1 a = +1, 



<7 = -l, 



This maximum at the stationary point can be seen in 
FigJSfl 

If, rather than a single electron, a beam of electrons 
is sent towards the line defect, the scattered electrons 
would be both valley- and spin-polarized. The former 
valley polarization, defined as 



V T 



(19) 



is shown in Fig.|4j This polarization is very close to 
V v = s'ma predicted previously in the absence of the 
local magnetic moments^ Similarly, we can define the 
spin polarization of the beam of transmitted electrons, 
as 



J2t T t,1 - J2r T t,- 



(20) 
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FIG. 4. Spin and valley polarization of the transmitted por- 
tion of a beam of electrons at the Fermi level after scattering 
off the line defect at the angle of incidence a. 



Although there is some spin polarization, as can be seen 
in Fig. 21 this polarization is small compared to the valley 
polarization. 



IV. CONCLUSIONS 

Graphcnc is a promising material for controlling elec- 
tron motion through scattering off well-defined defects. 
Herein, we have shown that electrons near the Fermi level 
scattered off an observed extended line defect can be both 
valley- and spin-polarized. The spin polarization, arising 
from local magnetic moments on sites adjacent to the line 
defect, is found to be less than 20%. The valley polar- 
ization, on the other hand, can reach near 100%. 

The valley filtering taking place at the line defect, 
which is very different from other suggested methods 
for obtaining valley polarization^ - — can be understood 
from its reflection symmetry. Consider an asymptotic 
graphene state near the Fermi level, which could be ex- 
pressed as |$ T ) = (| A) +ie^ l6 \B)) /y/2, where \A) and 
\B) refer to the two graphene sublattices and 9 is a pseu- 



dospin angle providing the phase relationship between 
the two sublattices. The only asymptotic graphene states 
also eigenstates of the reflection operator, which maps the 
A sublattice on one side of the line defect onto the B sub- 
lattice on the opposite side, and vice versa, are those with 
6 = ±7r/2. These states are symmetric and antisymmet- 
ric, respectively. Antisymmetric states must have a node 
at the line defect, making the coupling across the line 
defect, and concomitantly the transmission, small. Al- 
though, the potential describing the local magnetic mo- 
ments is a source of scattering, the transmission through 
the symmetric states is generally good. As a result, the 
line defect allows Bloch waves with 9 = tt/2 to transmit, 
while blocking those with 9 = —tt/2. 

The asymptotic graphene state |$ T ) is, in general, 
not an eigenstate of the reflection operator and has 
9 =/= ±7r/2. It can, however, always be written as a su- 
perposition of the symmetric and antisymmetric states. 
If the transmission through the symmetric and antisym- 
metric states are 1 and 0, respectively, the transmission 
probability can be obtained from the modulus square of 
the symmetric component of the incident graphene state, 
i.e. (1 + sin#)/2. From the graphene eigenstates, it can 
also be shown that the pseudospin angle 9 = roe, yield- 
ing a transmission probability equal to that in Eq. (|16|) in 
the absence of the potential describing the local magnetic 
moments £ 

To summarize, the valley filtering is a consequence of 
the imbalance between the transmission probabilities for 
the symmetric and antisymmetric components of the in- 
cident graphene state. This imbalance originates from 
the symmetry of the line defect structure. As neither the 
introduction of longer-range interactions,-^ distortion^ 
or potentials from the presence of local magnetic mo- 
ments, considered herein, can undo the imbalance, we 
conclude that the valley filtering is a robust property of 
the graphene line defect for high angles of incidence. 
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